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Abstract 

We use the spin-rotation-invariant Green’s function method as well as the high-temperature 
expansion to discuss the thermodynamic properties of the frustrated spin-S J 1 -J 2 Heisen¬ 
berg magnet on the body-centered cubic lattice. We consider ferromagnetic nearest-neighbor 
bonds Ji < 0 and antiferromagnetic next-nearest-neighbor bonds J 2 > 0 and arbitrary spin 
S. We find that the transition point between the ferromagnetic ground state and the 
antiferromagnetic one is nearly independent of the spin S, i.e., it is very close to the classical 
transition point jp 1 ' 1 * = At finite temperatures we focus on the parameter regime 

J 2 < with a ferromagnetic ground-state. We calculate the Curie temperature Tq(S,J 2 ) 
and derive an empirical formula describing the influence of the frustration parameter J 2 and 
spin S on Tc- We find that the Curie temperature monotonically decreases with increasing 
frustration J 2 , where very close to the Tc (Wpcurve exhibits a fast decay which is well 

described by a logarithmic term l/log(||Ji| — J 2 ). To characterize the magnetic ordering 
below and above Tc, we calculate the spin-spin correlation functions (SoSr), the sponta¬ 
neous magnetization, the uniform static susceptibility xo as well as the correlation length £. 
Moreover, we discuss the specific heat Cy and the temperature dependence of the excitation 
spectrum. As approaching the transition point some unusual features were found, such 
as negative spin-spin correlations at temperatures above Tc even though the ground state is 
ferromagnetic or an increase of the spin stiffness with growing temperature. 


1 Introduction 


Heisenberg models with competition between nearest-neighbor (NN) bonds J\ and next-nearest- 
neighbor (NNN) bonds J2 (J1-J2 models) can serve as canonical models to study frustration effects 
in magnetic systems. In case of ferromagnetic NN exchange, Ji < 0, and antiferromagnetic NNN 
exchange, J 2 > 0, the ferromagnetic ground state gives way for an antiferromagnetic state with 
zero net magnetization at a critical value that depends on the dimension D of the system m- 
In lower dimension D < 3 and for low spin quantum number S the ground state for J 2 > ,/f may 
have unconventional properties, see e.g. Refs. pUffllElE]- Although, in the ferromagnetic regime, 
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J 2 < J|, the ground state is the simple fully polarized ferromagnetic state, the influence of the 
frustrating coupling J 2 on the thermodynamic properties can be strong, as it has been recently 
discussed for systems with dimension D — 1 and D = 2 [SHH]- 

In the present study we extend our previous studies of the frustrated Ji-J 2 Heisenberg ferro- 
magnet in dimensions D — 1 (Refs. [HIl'Ss 3]) and D = 2 (Ref. [XT] ) to the three-dimensional 
case. In D = 3 the body-centered cubic (BCC) lattice is appropriate to compare with D = 1 and 
D = 2, since the antiferromagnetic J 2 bonds are in competition with Ji, but do not frustrate each 
other [T5Tn~8llX9] . 

The Hamiltonian of the spin-S' J 1 -J 2 Heisenberg model on the BCC lattice considered in this 
paper is given by 

H = J 1 Y,Si-S j + J 2 Y,S i -S j7 (1) 

(m) M 

where Sf — S(S + 1), and (i,j) runs over NN bonds and [i,j\ over all NNN bonds. We consider 
J\ < 0 and J 2 > 0. In what follows we define the energy scale by setting .J\ — —1. In the classical 
limit S —> 00 , the ground state is ferromagnetic for J 2 < J 2 ’ clab = ||Ji|, and it is a collinear 
antiferromagnetic state with the magnetic wave-vector Q = (ir, it, n) for J 2 > J 2 ’ clab [Tol ITS] . In 
what follows we want to discuss the thermodynamics of the model (P) for arbitrary spin quantum 
number S and focus on the parameter regime J 2 < J|, where the ferromagnetic ground state is 
realized. 

As in our previous papers on frustrated ferromagnets, mmm we use the rotationally invari¬ 
ant Green’s function method (RGM) to study the influence of the frustrating NNN bond J 2 and 
of the spin quantum number S on the thermodynamic properties of the model p. This Green’s 
function approach has been applied successfully to several quantum spin systems dnunumM. 
We complement the Green’s function study by using the high-temperature expansion (HTE) tech¬ 
nique [40l [4THT3ll44] . 

The paper is organized as follows: In Sec. 0 the RGM is developed for the model P and a 
brief instruction of the HTE approach is given. Some relevant features of the zero-temperature 
transition between the ferromagnetic ground-state and the collinear antiferromagnetic state are 
discussed in Sec. [3] Finite-temperature properties are then analyzed in Sec. 0 We summarize our 
findings in Sec. [5] 

2 Methods 

2.1 Rot at ion-invariant Green’s function method (RGM) 

As already mentioned above, the RGM is a well established method in the held of magnetic 
systems including strongly frustrated quantum spin systems. It was initially introduced by Hondo 
and Yamaji for spin-1/2 systems [T9] . Later on, an extension to arbitrary spin S was developed, 
see, e.g., Refs. [I X li2D lf221[Xfl S6] . 

To determine the spin-spin correlation functions and the thermodynamic quantities, we calcu¬ 
late the transverse dynamic spin susceptibility Xq“(o;) = — ((S'q ; S'Z q )) a , (here, ((...;.. .)) u denotes 
the two-time commutator Green’s function [X5HXT] ). Taking the equation of motion up to the sec¬ 
ond step (the decoupling of the higher order Green’s function is performed in the second-order 
equation of motion, i.e., one order beyond standard random-phase approximation) and supposing 
spin rotational symmetry, i.e., (Sf) = 0, we obtain cu 2 ((S' q ; S'Z q )) aj = M q + ((—5'+; <ST q )) w with 
M q = ( [[S q , H], ST q ] ) and — S+ = [[S q , H], H~\. For the model (P the moment M q is given by 
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the exact expression 


Mq — — 2 JiZiCooi(l — 7q^) — 2J2Z2Cl0l(l “ 7q^)j 
where c hkl = c R = {S£Sr) = 2(S 0 S fi )/3, R = ha i + ka 2 + la 3 , a t = \ V . [(1 - 26 e. 


( 2 ) 


[ej are 


the Cartesian unit vectors), 7 q = cos ^ cos ^ cos and 7q = ^(cosg x + cos q y + cosg 2 ). The 
quantities z\ = 8 and z 2 = 6 are the NN and NNN coordination numbers of the BCC lattice, 
respectively. For the second derivative — S+ we use the standard decoupling of the RGM, see, e.g., 
Refs. [f0lfII |]l^fl9l - l27l[35l - IF?] . To be specific, in —S) + we apply the decoupling in real space 


S?S+SZ = a iik (S?S,;)Sf + a Stk (S+Si)S? 
s+s;s+ = (s;st)s+ + X u(S+Sr)St, 


(3) 

(4) 


where i 7 ^ j 7 ^ k 7 ^ i. The quantities (x li3 and X t .j are vertex parameters introduced to improve the 
decoupling approximation. Note that products of three spin operators with two coinciding sites, 
see Eq. (QJ, appear only for S > 1, i.e., the vertex parameter X ltJ is irrelevant for S — 1/2 [T3ll22U36j . 
By using the operator identity Sf = S+S~ — Sf + (S)) 2 the expectation value (S',”S') 1 ") entering 
Eq. (J4]) is given by 

(S"5+> = (S+S-) = |S(S + 1), (5) 


where (S')) = 0 was used. For systems with a ferromagnetic ground state it is known from 
Refs. [Illl30ll35l - [37] that we may set a i)k = a and A it j = A, i.e., we approximate the various 
multispin correlation functions appearing in the equation of motion on an equal footing. We 
obtain — S+ — WqS'q and 


Xq M 


-((S+; SZq)) 


Mg 

^q “ ^ 


2 ’ 


( 6 ) 


with the spin-wave excitation dispersion relation 


Wq = ^l(7q ) -l) Q i (c001-C002 -3 (ci 0 1 + C 121 )+^lC 00 l7^ ) ) 

— 1) + ^^(Xq'’ — 1)) ^ c 001 + gS'(S' + 1)^ 

+ J2^2(7q 2 ^ — l) a ( C 101 — c 202 — 4Ci21 + ^ClOlTq 2 "*) 

+JiJ 2 ^2(7q 2) - !)« (-4c 120 — 12cooi + ^iCioi7q } ) 

+ Jl J 2 ^l( 7 q 1 ' ) — l)a (— 3 ci 20 — 3Cooi + ^2C00l7q 2 ' ) ) 

+ J1J2Z1Z2C00I { 7 q^ — l) (hq) — l) ) (V 


where we have used, on grounds of symmetry, coi_i = C 121 and cn_i = C 221 = < 7.20 • 

The correlation functions Cr = jj Cqe iqR are determined by the spectral theorem, 


Mr, 


Cq = (S+^_ q ) = ^-[1+ 2n(Uq)}, 


( 8 ) 


where n(cu) = (e u ^ T — 1) 1 is the Bose-Einstein function. Taking the on-site correlator cr = 0 and 
using the operator identity ([ 5 ]), we get the sum rule 


^£s, = fs(s + i). 

q 


(9) 
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To calculate the uniform static spin susceptibility xo — Xq=o(^ = 0), we use the relation x q (cu) = 
Xq Z (u;) = ly+ _ (u;). In the limit u = 0 and q = 0 we get from Eq. (J6]) 


where 


Xo = 


Cooi^l + C 101-/2 


A 


XO 


( 10 ) 


A xo = ^i a ( 9 c ooi — C002 — 3 (cioi + C121)) — </i</2«(cooi — 8 cioi + 7 ci 2 o) + ^l a ( 7 c 101 — 4 ci 2 i — C220) 

— A (J1C001 + ^ 2 c on) — 2 *S'(>S' + 1) (J/ + Jf) . (11) 

The phase with magnetic long-range order at T ^ Tc is described by the divergence of the static 
susceptibility at the ferromagnetic ordering vector Q 0 = 0, i.e., by Xo 1 — 0. In the long-range 
ordered phase the correlation function cr is written as [2T[l2H 361139] 

Cr = v E c i e “' R + c ’ < 12 ) 

q(^o) 


where c q is given by Eq. (T51) . The condensation term C determines the magnetization M which is 
given in the spin-rotationally invariant form by 


M 2 




(13) 


According to Eq. ca the condition for ferromagnetic long-range order is finally given by A Xo = 0, 
where A Xo is dehned in Eq. dill) . The magnetic correlation length above Tc is obtained by 
expanding y q in the neighborhood of the vector Qo [THESES]- For the ferromagnet we have 
Qo = 0 , and we find 


£ 2 

A 


A + B 16a (cqq-lJ) 5 + A 01 
16 (cooi Ji + C 101 J 2 ) 7A X0 

— J 2 J2<y ^12Cq 01 + Cqoi(31Cioi — 4ci2o) + Cioi(cqo 2 + 3(cioi + Ci2i))^ 


B 


-001 


— C101 f Aq 
= —3JiJ%a 
+ Cooi ( Aeon 


+ j5(5+l)V 
cooi(36cioi — 4ci2i 
+ §5(5+l)V 


C 220 ) + cioi(16cioi + 3c 



(14) 


(15) 


(16) 


To evaluate the thermodynamic properties, we have to determine the correlation functions cr, 
the vertex parameters a and A, and, for T < Tc, the condensation term C. For the correlation 
functions we use Eq. (fT2]i in combination with Eq. (j8jl . Further we use the sum rule (J9]) and the 
condition for long-range order, A Xo = 0, relevant for T < T c . That is sufficient for S' = 1/2, where 
the vertex parameter A is not present. For S > 1/2 we need one more equation. For T = 0 we use 
the known correlation functions of the ferromagnetic ground state, Cr{T = 0) = | S5r 0 + | S 2 , to 
obtain a(T = 0) = 3/2 and A(T = 0) = 2 — 1/S. To find the missing equation for S' > 1/2 and 
T > 0, we follow Refs. [3111551136] and consider the ratio 


r(T) 


A(T) - Aco 

a(T) 


r(0) 


(17) 
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as independent of temperature. Here A^ and are the vertex parameters for T —> oo which 
can be easily determined as = 1 — 3 [45(5 + l)]” 1 and a^ = 1. To solve the set of RGM 
equations numerically we use Broyden’s method, [48] which yields the solutions with a relative 
error of about 10” 8 on the average. The momentum integrals are done by Gaussian integration. 

In the low-temperature limit, several analytical expressions can be found from the RGM equa¬ 
tions which are identical to those obtained by the linear spin-wave theory (cf., e.g., Ref. [47]): 
The dispersion relation at zero temperature and |q| -C 1 (long-wavelength limit) is given by 
cu q (T = 0) = p s (0)q 2 , where p s (0) = p s (T = 0) = 5|Ji + J 2 I is the spin stiffness. Note that, 
by contrast to the one- and two-dimensional frustrated J1-J2 Heisenberg ferromagnet, 10. 11,13] 
the stiffness remains finite at the classical transition point = ||Ji|. The spin stiffness also 
determines the low-temperature regime of the magnetization and the specific heat, 


M 

~S 


c 


v 


c(i) 

( T V 

(18) 

25 

\47rps(0)/ 

15C(S) / 

f T \ % 

(19) 

8 ' 

v 4 Wb(0)/ 


where <C(s) denotes the Riemann zeta function. 

Some remarks are in order here with respect to the comparison between the RGM and the 
standard random-phase approximation (RPA), see, e.g., Refs. [35^ 4511491151] . The spin-wave ex¬ 
citation energies calculated within the RGM, see Eq. ([7]), exhibit a temperature renormalization 
that is wavelength dependent and proportional to the correlation functions, so that the existence 
of spin-wave excitations does not imply a finite magnetization. That is in contrast to the RPA, 
where the temperature renormalization of the spin-wave excitations is independent of wavelength 
and proportional to the magnetization, see, e.g., Refs. [45][46]. Correspondingly, in our theory 
the stiffness becomes temperature dependent p s (T) = + J 2 ) (cooi^/i + C 101 J 2 ) due to the 

temperature dependences of cc, C 001 , and cioi- Moreover, it is known [35 1 145 1 146] that the RPA fails 
in describing magnetic excitations and magnetic short-range order for T > Tc, reflected, e.g., in 
the specific heat. On the other hand, the magnetization for T < Tc is typically well described 
within RPA. 


2.2 High-temperature series expansion 

A very general approach to calculate thermodynamic quantities at high and moderate temper¬ 
atures is the high-temperature expansion (HTE) technique [401444] . Here we use a recently de¬ 
veloped tool |S| to calculate the HTE series of the specific heat and the uniform static suscep¬ 
tibility of model (Op) up to order ten. The region of validity of the HTE series can be extended 
by Pade approximants [401 - 1441153] . The Pade approximants [m, n] are ratios of two polynomials 
[m,n\ = P m (x)/R n (x ) of degree m and n and they provide an analytic continuation of a function 
g(x) given by a power series. As a rule, Pade approximants with m ~ n provide the best results. 

We can use the HTE series for the susceptibility Xo = Yh n c nT~ n to calculate the Curie tem¬ 
perature. Here we use two variants to extract Tc from the HTE series. 

(i) We analyze the quotient q n = c n /c n -i, see, e.g., Refs. [441154H56] . For a three-dimensional 
ferromagnet the critical behavior of Xo is given by Xo(T) oc (T — Tc)~T The expansion of 
(T — Tc) -7 in powers of 1/T yields for the quotient q n = Tc + (7 — l)yy- For higher orders n the 
HTE series of the model (]!]) obeys this relation, supposed that a finite critical temperature exists. 
We made a fit of our HTE data for q n linearly in 1/n including data points for n = 5,.. ., 10 to 
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get approximate values for Tq and the critical index 7. The mean square deviation of the linear 
fit provides information on the reliability of this estimate. 

(ii) A more sophisticated and more powerful approach to determine Tc and 7 is based on the 
so-called differential approximants (DA) [5TTf6Tj . Within this method one considers differential 
equations of the form Ylt=o Sv{P)Xo\P) + D(/3) = 0, where S v ((3) and Y((3) are polynomials in 
(3=1 /T and xb^ denotes the z/-th derivative of the susceptibility xo- (Note that the usual Pade 
approximants follow from K = 0.) The polynomials S„((3) and Y((3) can be chosen arbitrarily with 
the constraint that the total number of free coefficients of all polynomials is equal to the highest 
order n of the HTE series (here n = 10). The coefficients of the polynomials can be determined by 
inserting the derivatives of the HTE series of Xo hr the equation given above which yields a simple 
linear system of equations for the coefficients. Assuming critical behavior for the susceptibility Xo > 
we can determine the Curie temperature Tc = 1/(3c from Sk{(3c) — 0 and the critical exponent 7 
is given by 7 = K — 1 — ■ By varying the degree of the polynomials S u ((3) and Y (/3 ) a set 

of 21 different differential equations was obtained which yields a corresponding set of values for 
Tc and 7. However, some of the differential equations lead to unphysical results (negative or even 
complex solutions for Tc) and have to be discarded. The average of the data for Tc and 7 obtained 
from the remaining differential equations yields typically precise estimates of Tc and 7, where the 
mean square deviation provides a measure of the accuracy (for more details see Refs. pTTI ETT] ). It 
is in order to mention here that the HTE naturally works well only if the Curie temperature is 
not too small. In particular, for increasing frustration J 2 towards the critical value one may 
expect that the critical temperature becomes small, and then the determination of T c and 7 from 
the HTE series fails. 

We mention that in Ref. [16] the HTE series for the staggered susceptibility was used to 
determine the Neel temperature Ty as a function of J 2 /Ji for the model ((TJ) with antiferromagnetic 
J\ . Moreover, the Curie temperature of the unfrustrated BCC ferromagnet (i.e., J 2 = 0) for spin 
quantum numbers S = 1/2,1 and 3/2 was calculated from the HTE series in Ref. [62] . 


3 Zero-temperature properties 

First we use the RGM formalism to study the zero-temperature properties. For J 2 < the ground 
state is the simple ferromagnetic eigenstate with (S'o-S'r) = S 2 (R/0). The excitation spectrum 
(u q (T = 0) depends on J 2 , see Fig. [0 Since the spin stiffness is given by p s (0) = S\ J\ + J 2 |, 
see Sec. 12.11 the excitation energy in the long-wavelength limit is reduced by frustration, but the 
leading q 2 -term in (u q remains finite even at J 2 = J|. In the classical model the transition to 
the collinear antiferromagnetic ground state with magnetic wave vector Q = ( 7 r, 7 r, 7r) is related 
by the emergence of a soft mode at q = Q, i.e., ojq = 0 at J 2 = J 2 ,clas , for ferro- as well as 
for antiferromagnetic J\. We mention here the difference to the Ji-J 2 ferromagnet in dimension 
D = 1 (Refs. I lTi . IXFll ) and D = 2 (Ref. [H]). In D = 1 there is no soft mode, the spin stiffness 
p s vanishes at J 2 = = |Ji|/4 (i.e., the leading term in (u q is proportional to q 4 at J 2 = J 2 ), 
the transition point is independent of S, and the magnetic wave vector Q changes smoothly 
from Q = 0 to Q 7 ^ 0 when crossing the transition point Jf, see also Refs. . In D = 2 the 

behavior is again different. At the classical transition point J 2 ,chlb = |-/i |/2 all excitation energies 
o; q with q = (g, 0) and q = (0,g), 0 < q < n, become soft (i.e., there is a flat part in o; q ). 
However, this flat soft part is irrelevant for the quantum model, since J((' quant < J^ ,cas i i.e., the 
transition between the ferromagnetic and the antiferromagnetic ground states is shifted to smaller 
frustration strength (j^’ quant ~ 0.41 J x \ for S = 1 / 2 , see Refs. P,G5]). 
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r p h n r 

Figure 1: Dispersion relation of the spin-waves at T = 0 and for different values of J 2 and fixed 
Ji = —1. The points in the Brillouin zone are defined as T = (0,0,0),P = (7r,7r, 7t),H = 
(0, 0, 2tt), N = ( 7 r, 7T, 0). Note that for T = 0 the dispersion w q / S is independent of S. 

In view of these results for D — 1 and D = 2 the question arises, whether for D = 3 the 
transition point of the quantum model is different from that of the classical model. Let us first 
mention that for the model with antiferromagnetic NN bond, J\ >0, there is a slight shift of to 
~ 0.7Ji > J 2 ,clas for the S' = 1/2 model jTol 1X71149] . To find for our model with ferromagnetic 
Ji we extend our RGM approach to the antiferromagnetic region, i.e., to J 2 > J%, see Fig. El 
We follow Ref. [5] and use the intersection of the ferromagnetic ground-state energy Efm/N = 
—41 ^ | -S' 2 + 3 J 2 S 2 and the RGM energy of the collinear antiferromagnetic state to estimate J|, 
see Fig. El For S = 1/2 we get ~ 0.63|Ji|. Since the RGM description of the collinear 
antiferromagnetic state with only one vertex parameter a is a poor approximation, [2TU231139] 
a more accurate value for can be obtained using an improved description of the collinear 
antiferromagnetic ground state with two vertex parameters an and 02 , see, e.g., Refs. [2Tll271l33ll39j . 
The RGM with two vertex parameters yields ~ 0.68| J\ \ for S' = 1/2 and reproduces the classical 
value J 2 ,clab = 2|Ji|/3 in the large-S 1 limit [63]. Note that the RPA value for J|, estimated by the 
vanishing of see the next section, is = 0.6799| Ji\ |6l], 

4 Finite-temperature properties 

4.1 Spin-spin correlation functions, magnetization and Curie temper¬ 
ature 

The spin-spin correlation functions (SoSr,), the magnetization M and the Curie temperature Tc 
represent the basic quantities to characterize the ferromagnetic phase. First we discuss the Curie 
temperature Tc as a function of the frustration parameter J 2 as shown in Fig. [3] for S = 1/2. 
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Figure 2: Various ground-state spin-spin correlation functions (SoSr) as a function of J 2 for the 
spin-1/2 model in the ferromagnetic (Q fm = 0) and the antiferromagnetic regime (Q afm = P = 
(7r, 7r, 7r)). The inset shows the intersection of the energies per site of the ferromagnetic state and 
the collinear antiferromagnetic phase that is used to estimate J|. For comparison the GS energy 
of the classical model is also shown (labeled by ’AFM cl’). 


(Note that HTE results for the Curie temperature for S — 1/2,1,..., 5/2 and selected values of J 2 
are also given in Table [2] in Sec. 14.21 1 We compare the RGM data with the FITE results, variants 
(i) and (ii), see Sec. 12.21 and the RPA results of Ref. [64], Moreover, for the unfrustrated BCC 
ferromagnet (J 2 = 0) we can compare with HTE data of Ref. {62] for S = 1/2, 1 and 3/2. Note 
that our HTE data for J 2 = 0 within the accuracy of drawing in Fig. [3] coincide with those of 
Ref. [62], where HTE data up to orders 14 (5 = 1/2), 12 (S' = 1), and 9 (S = 3/2) were used. For 
J 2 = 0 and S — 1/2 the RGM (RPA) value is about 7% (13%) larger than the HTE value. The 
general dependence on J 2 is very similar for all methods. As expected, the error bars of the HTE 
estimates of T ( ; become larger with increasing J 2 , and the HTE fails as approaching J|, since the 
relevant temperatures are too small. Obviously, the more powerful approach using DA (i.e., HTE 
variant (ii)) works more reliable until large values of J 2 . By contrast to the HTE approach, the 
RGM and RPA provide 7% data until the classical transition point J 2 ,clas . As approaching , 
the Curie temperature remains quite large lip to about J 2 = 0.65 and finally drops down quickly 
at the transition point. The fast decay very close to J 2 ’ cla& is well described by a logarithmic 
term l/log(||Ji| — J 2 ) (cf. also Refs. [36] and [65]). The inset in Fig. [3] demonstrates that the 
RGM data for the Curie temperature scaled by S(S + 1) are only weakly dependent on S. The 
S'-dependence of the HTE estimate of Tq/S(S + 1) is more pronounced, but it remains also weak 
(e.g., for J 2 = 0.3 the difference between the extreme quantum case S = 1/2 and the classical 
limit S —y oo is about 20%). Hence, the curves for S = 1/2 shown in the main panel of Fig. [3] 
may be more or less considered as representative for arbitrary S. A good description of the RGM 

























Figure 3: Main panel: Curie temperatures as a function of the frustration parameter J 2 

and J\ — — 1 for spin quantum number 5 — 1/2 obtained by the RGM, RPA (data from Ref. |64j). 
and by the HTE, variants (i) and (ii). Inset: Curie temperatures as a function of 1/5 for 

various frustration parameters J 2 . Note that within the RPA approach is independent of 

5. 

Table 1: Coefficients A and B of the empirical fit function of the Curie temperature Tc/S(S + 1), 
see Eq. (1201) . for various spin quantum numbers S. 



5=1/2 

5 = 1 

S = 3/2 

5 = 2 

5 = 5/2 

5 = 100 

A 

3.29 

3.42 

3.34 

3.27 

3.22 

3.04 

B 

1.41 

1.41 

1.35 

1.31 

1.28 

1.18 


data for T C /S(S + 1) in the entire range 0 < J 2 < || J\\ is provided by the fit function 


T c A(S) 

5(5 + 1) B(S) - log(|| A| - J 2 ) ’ 


( 20 ) 


The fit parameters A and B as a function S are given in Table [0 Obviously, the 5-dependence 
of A and B is weak. 

Next we consider the spin-spin correlation functions (SoSr,). In Fig. [4] we compare the tem¬ 
perature dependence of the correlation functions between nearest neighbors (R = (1,0,0)), next- 
nearest neighbors (R = (1,1,0)), and between spins separated by R = (2,2,0) for J 2 = 0 and 
J 2 = 0.6. For T < Tc the behavior of the correlation functions for J 2 = 0 and J 2 = 0.6 is quite 
similar, although their decay with increasing T is faster for larger S. Interestingly, above Tc the 
correlation functions for particular separations R for larger values of frustration become negative 
irrespective of the fact that the ground state is still ferromagnetic. These negative correlators are 
precursors of the collinear antiferomagnetic Neel order present for J 2 > J|, since they belong to 
separations R, e.g. R = (1,1, 0), connecting spins on different sublattices of the large-J 2 collinear 
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Figure 4: Spin-spin correlation functions (SoSr)/^ 2 as a function of the normalized temperature 
T/Tc for spin quantum numbers S' = 1/2 (red lines) and 5/2 (black lines) and for two values of 
the frustration strength J 2 = 0 (left) and J 2 = 0.6 (right). 

antiferromagnetic phase, cf. also Fig. [2j The general behavior of (SoSr) for T > Tc at large 

separation |R| is given by (S 0 Sr) ^ e~ (|R|~ D+2_ '' 7 ), where D = 3 and the critical exponent rj 
is close to zero, see, e.g., Refs. [66] [68]. The long-distance behavior of RGM correlation functions 
at T = Tq turns out to fit well to a power law with rj ~ 0 for arbitrary S and J 2 . 

In Fig. 0we show the temperature dependence of the magnetization for three values of S and 
for two values of the frustration parameter, J 2 = 0 and J 2 = 0.6. The influence of J 2 on the M/S 
vs. T/Tc curve is rather small, an increase of the spin quantum numbers S leads to a flattening 
of the M/S vs. T/Tc curve. 

According to Refs. [B5] and pni the spontaneous magnetization of a ferromagnet can be de¬ 
scribed by the formula 


M 

~S 



(1-0 



( 21 ) 


where 0 is the so-called shape parameter and j3 is the critical index of the order parameter, 
which is (3 = 1/2 within our RGM approach, cf., e.g., Ref. [36]. This formula yields the correct 
low-temperature behavior and a critical behavior for T —>■ Tc with the exponent /?. Indeed, a 
corresponding fit of our RGM data for M using the shape parameter 0 as a single fit parameter 
yields a very good description of M(T) in the entire temperature region 0 < T < Tc- We show the 
shape parameter (j) in the inset of Fig. [5] Based on the low-T properties known from spin-wave 
theories and on early estimates of T C) eh Kuz’min et al. [69],[7D] argued that the dependence of 
0 on S and J 2 should be weak. Moreover, they found a relation 0 < 0 < 5/2. Our results for 0 
support these statements. 
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Figure 5: Main panel: Magnetization M/S as a function of the normalized temperature T/Tc for 
spin quantum numbers S = 1/2,2 and 100 and for two values of the frustration strength J 2 = 0 
(solid lines) and J 2 = 0.6 (dashed lines). Inset: Shape parameter 0 of the fitting function (1211) in 
dependence on l/S for various values of frustration J 2 . 



1 1.5 2 2.5 3 

T/T c 

Figure 6: Inverse scaled susceptibility 1/(| J 2 + Ji\)x as a function of the normalized temperature 
T/T c . Main panel: RGM results for S = 1/2 (solid lines) and S = 5/2 (dashed lines) for 
frustration parameters J 2 = 0 (blue lines) and J 2 = 0.5 (red lines). Inset: Comparison of RGM 
data (lines) with Pade approximants [6,4] of the 10th order HTE data (symbols). 
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Table 2: Critical exponent 7 of the susceptibility and Curie temperature Tc/S(S ± 1) obtained 
from HTE (i) and HTE (ii), cf. Sec. 12.21 for different spins S and frustration parameters J 2 and 


J\ = -] 


7 

J 2 = 

0 

h = 

0.1 

J 2 = 

0.3 

S 

HTE (i) 

HTE (ii) 

HTE (i) 

HTE (ii) 

HTE (i) 

HTE (ii) 

1/2 

1.53 ±0.05 

1.34 ±0.04 

1.51 ±0.10 

1.33 ±0.08 

0.26 ±2.07 

1.38 ±0.40 

1 

1.49 ± 0.04 

1.44 ± 0.20 

1.42 ±0.10 

1.42 ±0.08 

0.46 ± 1.30 

1.39 ±0.10 

3/2 

1.46 ± 0.04 

1.40 ±0.10 

1.41 ±0.09 

1.41 ±0.07 

0.57 ± 1.18 

1.36 ±0.03 

2 

1.45 ± 0.04 

1.39 ±0.10 

1.39 ±0.09 

1.43 ±0.10 

0.61 ± 1.13 

1.40 ±0.10 

5/2 

1.44 ± 0.04 

1.45 ±0.20 

1.39 ±0.09 

1.40 ±0.07 

0.64 ± 1.09 

1.39 ±0.07 

Tc 

S(S+ 1) 

J 2 = 

0 

J 2 = 

0.1 

J 2 = 

0.3 

S 

HTE (i) 

HTE (ii) 

HTE (i) 

HTE (ii) 

HTE (i) 

HTE (ii) 

1/2 

1.67 ±0.01 

1.69 ±0.01 

1.50 ±0.02 

1.55 ±0.04 

1.41 ±0.42 

1.19 ±0.09 

1 

1.88 ±0.01 

1.88 ±0.03 

1.71 ±0.03 

1.70 ±0.01 

1.55 ±0.29 

1.33 ±0.01 

3/2 

1.96 ±0.01 

1.96 ±0.02 

1.79 ±0.03 

1.77 ±0.01 

1.60 ±0.27 

1.40 ±0.02 

2 

1.99 ±0.01 

1.99 ±0.02 

1.82 ±0.03 

1.80 ±0.02 

1.63 ±0.26 

1.43 ±0.01 

5/2 

2.01 ±0.01 

2.00 ±0.03 

1.84 ±0.03 

1.82 ±0.01 

1.64 ±0.26 

1.44 ±0.03 


4.2 Susceptibility and correlation length 

Next we discuss the uniform susceptibility and the correlation length £. In addition to the 
spin-spin correlation functions, see Sec. 14.11 both quantities characterize the paramagnetic phase 
above T c . We show the inverse susceptibility in Fig. [61 The susceptibility multiplied by (|./ 2 + ./11) 
shows an almost universal behavior as a function of the normalized temperature T/Tc , i.e., the 
dependence on S and J 2 is small. The RGM results agree well with the HTE data (see the inset 
in Fig. [6]). While the RGM critical exponents deviate from the correct values, the HTE approach 
provides more accurate exponents. We present HTE data for 7 in Table [2j As expected, 7 is almost 
independent of S and J 2 an d it is close to the correct value 7 ~ 1.39 of the three-dimensional 
Heisenberg ferromagnet, see, e.g., Ref. [68] . As already noticed in Sec. 14.11 the HTE approach (ii) 
to determine T c and 7 using differential approximants is more accurate and works well until large 
values of J 2 . Thus, the analysis of the quotient of the HTE series for xo- he. HTE (i), starts to 
fail already at J 2 ~ 0.3, see the large least-square deviations in Table [2] 

The RGM results for the correlation length are shown in Fig. [0 Again the influence of the 
spin quantum number S and the frustration J 2 on the £(T/T C )~ curve is weak. The decay of £ 
with temperature is fast, already at T ~ 1.37b the correlation length is of the order of the lattice 
constant. 
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Figure 7: Correlation length £ as a function of the normalized temperature T/Tc for S = 1/2 
(solid lines) and S = 5/2 (dashed lines) for frustration parameters J 2 = 0 (blue lines) and J 2 = 0.5 
(red lines). 



Figure 8: Temperature dependence of the dispersion relation cn q /S of the spin excitations for 
different values of J 2 and fixed J\ = —1. The points in the Brillouin zone are defined as T = 
(0, 0, 0). P = ( 7T, 7T, 7r), H = (0, 0, 27 r), N = (7T, 7T, 0). 
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4.3 Excitation spectrum and specific heat 

Let us now consider energetic quantities. By contrast to the RPA, the RGM provides a reasonable 
description of the excitation spectrum cn q for arbitrary temperatures. We show the dispersion 
relation for various temperatures, spin quantum numbers and frustration parameters in Fig. HI 
First we notice that for finite temperatures (u q / S becomes dependent on S, see Figs. [Hh and b. ft 
is also evident that for a wide range of temperatures up to about 50% of the Curie temperature 
the changes in cn q are quite small, see also Fig. |9j The main aspects of the influence of the 
temperature on cu q are the following: For small values of J 2 < 0.3 the excitation energies for all 
q-vectors are reduced by increasing T. Around the maximum at H = (0, 0, 2tt) this reduction is 
present for all values of J 2 . However, the temperature dependence of cn q around the dip/minimum 
at the soft-mode wave vector q = P = (tt, 7t, 7 r) as well as for q — y 0 differs for smaller J 2 from 
that for larger J 2 . (Note that the second dip/minimum between the H- and N-points at about 
q = 0.66(7t, 7 r, tt) is related to that at q = P, i.e. it appears only due to the choice of path H — * N 
in the Brillouin zone.) At q = P the excitation energy cup decreases with growing T for smaller 
J 2 , but it increases for larger J 2 , see Fig. [HE The influence of T on u; q near the T-point, i.e., at 
small |q|, is more subtle: For T < Tq the temperature dependence of the excitation energies is 
given by the temperature dependence of the spin stiffness p s describing the quadratic term in cu q 
for small |q|, while for T > Tc a linear term in tu q determines the behavior near the T-point. We 
show p s (T) for S = 1/2 for various values of J 2 in Fig. [9E- As expected for ferromagnets, [72l 174] 
the stiffness becomes smaller with increasing of T for J 2 < 0.6. Interestingly, for J 2 > 0.6 p s (T) is 
growing with T. This unusual behavior of p s (T) near the zero-temperature transition point to an 
antiferromagnetic state has recently been discussed in Ref. m using non-linear and self-consistent 
spin-wave theories. Note, however, that in Ref. ra the frustrated FCC ferromagnet is considered, 
where p s (T = 0) =0 at the transition point to the antiferromagnetic ground state, whereas for 
our BCC model p s {T = 0) > 0 at J? l . For small and moderate frustration, the spin-wave approach 
of Ref. [73] yields a decrease of p s with T according to p s (T) = p s (0) — BT % p = 5/2 and B > 0. 
As approaching the zero-temperature transition point to the antiferromagnetic ground state, the 
exponent p changes to p = 5/4 and the prefactor B becomes negative, i.e. the stiffness can grow 
with increasing T. Our RGM data for p s (T ) confirm these predictions: We find that p s (T ) is well 
described by the power law given above with B > 0 and p ~ 2.5 for J 2 < 0.5, whereas for J 2 = 0.6 
we have B & —0.05 and p ~ 1.5. 
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Figure 9: Temperature dependence of the spin stiffness p s (left) and excitation energy c op at the 
soft-mode q - vector P = (tt, tt, tt) (right) for spin S = 1/2. 



0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 


T/T c 


Figure 10: Specific heat Cy as a function of the normalized temperature T/Tq- Main panel: 
RGM data for S = 1/2,1, 3/2 and 5/2 for frustration parameters J 2 = 0 (solid lines) and J 2 = 0.5 
(dashed lines). Inset: Comparison of RGM data (dotted lines) with Pade approximants [6,4] of 
the 10th order HTE data (symbols). 

The specific heat Cy is shown in Fig. [TQ1 By contrast to the susceptibility and the correlation 
length the influence of the spin S and the frustration J 2 is noticeable. The Cy(T) curves show the 
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characteristic cusplike behavior at the transition temperature around Tq for small spin quantum 
numbers S indicating the second-order phase transition. With increasing S in the long-range 
ordered phase the slope of the Cy curves near Tq decreases, and the cusplike shape develops to 
a kinklike one. This behavior for larger S may be considered as a shortcoming of the RGM, cf. 
Ref. |36j . For T > T c we can compare the RGM data with HTE results, see the inset in Fig. (TUJ 
As for the susceptibility, the agreement is very good. The low-temperature behavior is given by 
Eq. ([T9]) , i.e., Cy(T ) increases with increasing J 2 . Note, however, that the specific heat as a 
function of the normalized temperature T/Tc shows the opposite trend, see Fig. fTOl 

5 Summary 

In this paper we use the rotation-invariant Green’s function method (RGM) to calculate thermo¬ 
dynamic quantities, such as the Curie temperature T c , the spontaneous magnetization M, the 
spin-spin correlation functions (SoSr), the uniform susceptibility y 0 ; the correlation length £, the 
specific heat Cy and the spin stiffness p s of a frustrated spin-S' Heisenberg magnet on the BCC 
lattice with ferromagnetic NN exchange J\ and antiferromagnetic NNN exchange J 2 . We focus 
on the ferromagnetic regime, i.e., J 2 < 2| Ji|/3, such that the ground state is ferromagnetic. For 
T > Tq our Green’s function approach is complemented by 10th order high-temperature expansion 
for the susceptibility and the specihc heat. 

The RGM goes one step beyond the random-phase approximation (RPA). As a result, several 
shortcomings of the RPA, see, e.g., Refs. [ TSlITSlITUHSTj . such as the artificial equality of the 
critical temperatures Tjv = Tc for ferro- and antiferromagnets or the failure in describing the 
paramagnetic phase at T > Tc , can be overcome. For the Curie temperature and the spontaneous 
magnetization M we derive simple fit formulas describing Tc as a function of S and J 2 and the 
temperature dependence of M(S,J 2 ). As approaching the ground-state transition point to the 
antiferromagnetic phase at ~ 2| JjI/ 3, the thermodynamic properties deviate from the ordinary 
ferromagnetic behavior. Thus, the spin-spin correlators may become negative at T > Tc indicating 
the collinear antiferromagnetic Neel order present for J 2 > J|, and the temperature profile of the 
spin stiffness p s for T < Tc exhibits an increase with T instead of the ordinary decrease. 

The present investigations are focussed on theoretical aspects, in particular, with respect to 
previous discussions of one- and two-dimensional frustrated ferromagnets. There might be some 
relevance for ferromagnetic compounds [75H77] especially near a quantum phase transition, e.g., 
driven by doping. 
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